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Detailed calculations of the contribution from off-shell effects to the quasiclassical tunneling of 
fusing particles are provided. It is shown that these effects change the Gamow rates of certain 
nuclear reactions in dense plasma by several orders of magnitude. 



Methods of field theory do not straightfor- 
wardly apply to systems at finite temperature. 
Naive perturbation theory fails in the infrared 
region of energy at high temperature but can 
be fixed by the partial HTL resummation [|TJ. 
This approach was generalized in |2| , and it was 
shown that fields without dispersion relation || 
are necessary for consistent formulation of the 
theory. 

The particle distribution over momenta, /(p), 
acquires non-Maxwellian tail in thermodynamic 
equilibrium [Q due to quantum effects, while the 
particle distribution over energies, f(co), remains 
Maxwellian. An interesting suggestion was made 
in Ref. ||, that the rates of various processes, 
including nuclear reactions, may be significantly 
increased by the quantum tail of the momentum 
distribution. Non-linear integral equation for the 
rate of inelastic processes as well as a final for- 
mula for the ionization rate was obtained in [Q. 
The formula for the nuclear reaction rates was 
found in || under assumption, that it is correct 
to substitute the particle momentum in the form 
e p = p 2 /2m into the known quasiclassical cross- 
section o~(uj) and then average it over the parti- 
cle distribution over momenta /(p) rather than 
the distribution over energies f{uS). This proce- 
dure becomes even more unclear , if we take into 
account the fact, that the momentum, p, and 
energy, u>, of a colliding particle are independent 
variables not connected by the usual dispersion 
relation 5{u - p 2 /2m) 

Note, that the averaging 
the distribution over energies 
Maxwellian, would give the 
rates || rather than the rates found in || , which 
are accelerated by many orders of magnitude in 
certain regimes. One can argue that aver- 
aging over momentum may not be completely 
wrong, because the dependence of the scatter- 



of o-{uj) over 
f(uj), which is 
usual reaction 



ing amplitude on energy can be neglected in the 
gaseous approximation P,|IO||. However, the con- 
nection between the scattering amplitude and 
the rate of nuclear reactions is not clear. Neither 
is it obvious that this argument can straight- 
forwardly lead to correct answer in the case of 
essentially non-linear processes of nuclear reac- 
tions. It was shown that non-linearities can be 
very subtle, subject to various non-trivial cancel- 
lations []llj] . It is therefore important to find the 
nuclear reaction rates from the first principles. 

In this paper we rigorously solve the quasiclas- 
sical problem of tunneling through the potential 
barrier. Our result demonstrates that the can- 



cellations of the type found in JO] do not play 
a role in the case under consideration. We dis- 
cuss this question in more details in the future 
publication |T2| . 

The tunneling particles undergo simultane- 
ous collisions with other particles of the plasma 
maintained in thermodynamic equilibrium. The 
plasma is assumed to be fully ionized. We 
use the Green-function technique, introduced by 
Keldysh and Korenman |H|,|l4j], and do not rely 
on any ad hoc assumptions || about any aver- 
aging procedure. Our final result for the nuclear 
reaction rate |15[ agrees with that postulated in 
0, apart from the coefficient depending on the 
masses of colliding and reacting particles. 

A number of different off-shell processes in 
plasmas were investigated within a framework 
of thermal field theory. Emission rate of soft 
photons from hot plasmas was calculated in [16 |. 
Using full spectral functions allowed to eliminate 
unphysical divergences and extend the range of 
validity of the results obtained by HTL resum- 
mation by taking into account off-shell parti- 
cle propagation. Landau damping was shown 
to lead to anomalous relaxation of the parti- 
cle distribution function 101 . Off-shell effects 
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due to particle collisions were treated within 
Born approximation in |Tj| . The relaxation rate, 
calculated with the help of full spectral func- 
tions, is slower approximately by a factor of 2 as 
compared to the rates calculated by solving the 
Boltzman equation [18]. This relaxation rate is, 



however, an integral quantity, insensitive to the 
details of changes of the kinetic Green function 
due to inclusion of off-shell effects. 

In contrast to the relaxation rate, inherent 
non-linearity of the tunneling problem makes 
the average nuclear reaction rate very sensitive 
to the behaviour of this function at large mo- 
menta. This explaines why we should expect 
off-shell effects to lead to significant corrections 
to the nuclear reaction rates. Indeed, we find 
that certain reactions are accelerated by many 
orders of magnitude as compared to the calcu- 
lations, performed using spectral functions with 
zero width |T]J . We solve the tunneling problem, 
neglecting Salpeter corrections due to screening, 
and assume that, the nuclear transformations are 
unaffected by the medium. 

Kinetic properties of the particles are de- 
scribed by the Green function G~ + (X 1 ,X 2 ), 
where X = (t,x) fI3 |. Equations for G h in- 



clude the interaction between the particles as 
well as their interaction with the external field. 
Therefore, the rate of change of this function, 
evaluated at any point after the barrier, will give 
the tunneling rate K(x, t) for the particles col- 
liding simultaneously with other particles |15| : 



K(x,t)= (d h -d t2 )G-+(t 1 x 1 ;t 2 x 2 ) 



X\^X2=X 



Now we will proceed with solution of the ki- 
netic equations for G~ + (Xi,X 2 ). They can be 
written as Dyson equations for G°[ 2 |13| : 



Ga/3 
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where we use subscripts 1..4 to denote depen- 
dence on X1..X4. Dependence on Gf 2 alS can be 



eliminated |L3[ by acting on both sides of the 
Eq. (0) with the operator 



cr 1 



.d_ Aj_ 

dti 2m 



where Af = + d 2 . Since we are interested 

J- 1 y 1 

in the steady state solution, one can see that, 
the properties of this operator can be fully ac- 
counted for if we introduce the function g(z, k), 
such that 



2m dz 2 



U{z) 



g(z,k) = -e k g(z,k), (3) 



where e k = k 2 /2m. Since the operator L is 
Hermitian, we will use the property of its eigen- 
functions g(z, k): 



g(z,k 1 )g*(z,k 2 )dz = 2nS(k 1 ~k 2 ) (4) 



Now let us solve for G^ + , which satisfies 



r 12 



0. 



(5) 



It is easier to find G^ 2 + by using its definition 



G (o)— 



r 12 



(6) 



where ip\ ip are the Heisenberg operators of cre- 
ation and annihilation of the particles at a point 
X, and < .. > means quantum and statistical av- 
eraging. Since ip evolves according to G~ l ip = 0, 
one can see, that it is equal to 

4>{X) = £ a^ k g(z, k)e-^ t+l ^ x± (7) 

q x fc 

Here a q±/ t is the annihilation operator of the 
particle with momentum (qj_,fc), and e q± fc = 
)q 2 ± /2m + k 2 /2m. 

We substitute ip(X) from (0) into Eq. (||) and 
obtain 



r°° dk 
2m / — 

J-oo 2ix 



5(u - e q±fc + /j l )n(q±k)g(zi, k)g*(z 2 , k), (8) 

where we denote (a qx fcflq ± fe) = n(q±k); /i is the 
chemical potential. 

It is easy to see that, in thermodynamic equi- 
librium we are considering, the function n(q_|_fc) 
is equal to Fermi distribution rip(e q± fc). To con- 
vince ourselves we take G^~ + (zi, z 2 ) at such 
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values of z\,z 2 , when U(z) is negligible and 
g(z, k) =>■ exp(ikz). One more Fourier transform 
Jexp[-iq(z 1 - z 2 )]d(z 1 - z 2 ) will give G^ + , 
which has to coincide with the known expres- 
sion = 2™% - e q + //)n F (e q ). 

Before proceeding further with solving for 
G- + (X U X 2 ) we first find K { ®\ which is the tun- 
neling rate obtained by using G^ 2 ~ + in Eq. (jl|). 
Since we are interested in behind the bar- 
rier, where U(z) = 0, we can exchange dt ± for 
iAi/2m as is clear from the form of the operator 
Gi 1 and Eq. @. Then the formula for the rate 
K( ' takes the form: 



v ' 2m J 2vr 7 



(2tt) 2 



Q{2 u wqj 



(*,0 



(9) 



where z 



2l - ^2- 



«l+22)/2, C 

12 through 
.. q v 2 ! C) with respect to 



In deriving @ we expressed G^ 0; 
its Fourier transform G>?2 



'fast" space rj_ = — r 2 and time t = t\ 



t 2 

as pre- 



variables and took the limit r, r 
scribed by Eq. ([I]). 

Finally we need to know the function g(z,k) 
in the region behind the barrier. We will use the 
standard quasiclassical expression |19f 



g(z, k) 



Gh 



q(z,ky/- 



; exp 



i / q(z', k)dz 



C k = A k exp 



z b (k) 



q(z', k)dz 



(10) 



(12) 



q(z,k) = J2m(e k -U(z)) 



The square of the exponential factor in Eq. ( p]) 
is the tunneling coefficient, W{k) |19| . We multi- 
ply it by a factor A k = (S(e k ) f^k) 1 , which takes 
into account nuclear physics effects, not consid- 



ered here. S(ek) is the astrophysical factor p0 

To find K^> from Eq. @ we substitute 
Eqs. ©-(Q) in Eq. ©, use z x = z + (/2, 
z 2 = z — (/2 and take d 2 derivative in the quasi- 
classical sense keeping only zero-th order terms. 
The result of the differentiation is 



a 2 



»72 [g{z 1 ,k)g*{z 2 ,k)] 



-q(z,k)C 2 k (13) 



JC=o 



Note, that q(z, k) = k is true behind the barrier. 

The problem of tunneling collision between 
particles m\ and m 2 can be reduced to the one- 
dimensional motion of a particle with reduced 
mass m r = m\m 2 j (mi + m 2 ) in the external po- 



tentianl U(r) ||19| . We make this reduction in 
the original Hamiltonian and then second quan- 
tize it. 

We need to perform calculations of G^~ + and 
K similar to those outlined above. However, 
now we have to do 3-D calculations, using spher- 
ical coordinates. So we replace L(,2)-operator 
with L x = A x /2m r , denote its radial part in 
spherical coordinates by L r and introduce the 
function g(x) according to 



#( x ) 



k 2 dk 
W) 



-,9kW 



9k(x) =J29ki{r)Yt 



tin 



I'm 



(L r -t/ { (r))^«(r) 



-^9ki[.r) } 



(14) 
(15) 

(16) 



where Yi m (9<j)) is a spherical harmonic and 
Ui(r) = U(r) + 2^5- Quite analogously to 
Eq. ([?D we express ip through g(x) as 



fclilTTll 



(17) 

Now using Eqs. ©, (|16|) , dl7|) and taking into 
account the contribution only from spherically 
symmetric modes, we obtain G^°^ + 



a (0) - 



1*1, *2 



2vri 



k 2 dk 
o (2vr) 3 



<5(u; - e k + ^)njr(e k )(/ fc (xi)flrjj(x2). (1* 
Rewrite the rate .ff , Eq. @, as 
i r du 



2m r J 2k 



a x1 g£°)- + ( Xi ,x 2 ) 



(19) 
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In the WKB approximation we get for the func- 
tion <7fc(r) the following formula 



9k(r) = -x(r) 
r 



X(r) 



q(r, k) 1 /' 



: exp 



i / q(r , k)dr' 

J To 



W(k) = exp 



q(r', k)dr' 



ro 



Q 



(r, k) = sj2m r (e k - U{r)), 



(20) 



(21) 



(22) 



(23) 



where the bar denotes dependence on the re- 
duced mass, e k = k 2 /2m r , C\ = A 2 k W(k), 
a = ZiZ 2 eie 2 and A\ = r^S^k)/^. We evalu- 
ate the derivative of Eq. fll9|) in the quasiclassical 
sense 



[^( X 2)A Xl ^fc(Xi)] 



xi^X2=ro 



r<2 



(24) 



and obtain the final result for 



r ^ri\ f 00 Ank 2 dk k S(e k )„ T/7 . 



(27T) 



m r e k 



(25) 



where nju(^) is the Maxwell distribution which 
depends on m r . 

It agrees with the answer obtained in || 
by averaging the quasiclassical tunneling fac- 
tor W(k) over the Maxwell distribution, 
(nM(v)a(m r v 2 /2)v). Hence, we see, that the dis- 
persion relation 8{uj — e q ) was implicitly used in 
||, which corresponds to an approximation of 
instantaneous, two-body collisions. Now we will 
not make this assumption and proceed with find- 
ing the tunneling rate, K, of the particles, which 
always collide and hence are never "in between 
collisions" . 

By using |I3|J7|, we find 



T^WH (G^{z) - Gi^zj) , (26) 



where we denote T, R — S A = «7 wq (z) = 
*7wqx(?>2)- Note, that 7^(2) is a non-linear 
function of G aP We now analyze Eq. © 

qualitatively, which will help us to find that part 
of the solution, which makes the largest contri- 
bution to K. 

In the region away from the barrier we can ne- 
glect the LHS of Eq. (p6|), and obtain: 



-n F {uj i 



R , 



(27) 



Therefore, Eq. (|26|) will allow us to propagate 
this solution into the region behind the barrier. 

We can find the first integral and formally inte- 
grate Eq. ( PE| ) along trajectories, which will lead 
to a sum of solutions of homogeneous and "in- 
homogeneous" equations (due to the last term 
in (p6|) ). As can be seen from Eq. (p6|) , "inho- 
mogeneous" solution will involve the integral of 
the product of j U q(z) and G^z) or G R * = G A . 
Qualitatively, it means, that in the region of 
z behind the barrier it will be proportional to 
a product of at least two exponential factors 
W{k). This is so, because both 7a, q (z) and 
G R (z) will depend on the integral of the prod- 
uct g(z±, k)g*(z 2 , k) and such product is oc W(k) 
as can be seen from Eqs. (]T0|) - (pT2|) . 

Explicitely, the expression for the retarded 
Green function can be obtained through the 
usual resummation IPI 



Gujqf^z) 



uj — e r 



(2i 



Since the advanced function, G A , is related to 
G R by G A = G R *, we find that their difference 
is 



G R — G A = 21m 



1 



7 



(G R y 2 + (imz R y 



(29) 



Now we will show that 7 is oc W(k), which will 
allow us to neglect the "inhomogeneous" term in 
Eq. (J26|), as being second order in W(k) since 
qr _ qA j g a j go p r0 p 0r tional to W(k). First, 
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note that Eq. (15) is written in the Furry rep- 
resentation. The Furry representation is needed, 
since we are interested in the behaviour of all the 
functions in the external field U(z). This means 
that the diagrammatic expansion of the opera- 
tor 7 = 2ImTt R in Eq. fl26|) is built upon zero-th 
order functions, which incorporate dynamics in 
the external potential and hence are exponen- 
tially suppressed in the region behind the bar- 
rier. This is in contrast to the usual expan- 
sion which utilizes zero-th order functions of free 
particles. (Such representation is usually called 
Furry representation in the literature.) In other 
words, these functions resolve the operator 




U'(z. 



cT 
dq j 



(30) 



which is the left hand side of Eq. fl26|) . There- 
fore, they are exponentially suppressed in the 
region of z behind the barrier. This can be 



also seen from the explicit expressions for the 
Green functions. The function G^~ + , is writ- 
ten in Eq. (§), and is proportional to the prod- 
uct g(zi, k)g*(z 2 , k). As discussed above, this 
product is exponentially suppressed behind the 
barrier. Analogously, the function G^ + ~ can be 
obtained from G^~ + by substitution n(qj_k) —> 
(1 — n(q±k)), and is also suppressed exponen- 
tially, since its dependence on z is the same. This 
is an important property, which we will make use 
of, while analyzing the behaviour of 7 in the re- 
gion behind the barrier. 

Now we can write explicitely expression of the 
first Born diagram for 7ojq x (^i, z 2 ) in the Furry 
representation from the Fig. 1 b) of Appendix A. 
There is a slight difference from the case of free 
zero-th order functions. Namely, there are addi- 
tional integrals over intermediate z 3 , ^-points in 
the loop (see Appendix A for notations): 



^pj^i^a) = / G a + (u,p ± - q ± ; z u z 2 )G£ (u, p±; z 4 , z 3 )G b p^ - q x ; z 3 , z 4 ) 



V(q±; z\ - z 4 )^(q_L; z 2 - z 3 )^^^-^^^dz 3 dz 4 , 



(2vr) 3 



(31) 



Since all the Green functions in Eq. (|3ll) are 
exponentially suppressed behind the barrier, we 
will obtain the largest contribution, if we re- 
strict integration over z 3 and £4 to the region 
away from the barrier, z > z*. The func- 
tions G£~ , G b + are not exponentially suppressed 
there. However, note that (zi, z 2 ) depends 
on zi, z 2 through the function G~ + in Eq. (|3T|). 



One can see that behaviour of £ (q, z) in the 
region of z behind the barrier will be determined 
by the behaviour of G~ + (u,pj_ — q_i_; 21,^2), 
Wigner-transformed as in Eq. (j32[). The same 
argument applies to which will be deter- 

mined by G+~ ■ Therefore, 7 = £~ + — 
which is 



= «*Hrf* " *)] (PZLM - ^M) - ») (32) 



will depend on the behaviour of these Green 
functions behind the barrier. 

But we know that, since both G^~ + and 
£(o)+- reso i ve the operator (p0|), they are expo- 
nentially suppressed behind the barrier. (This 
is also proven by their explicit expressions, as 
explained above). Therefore, 7 a;q± (g, z) is also 
exponentially suppressed there. 

More formally, we can obtain explicit expres- 



sion for the 7 as far as its dependence on z is 
concerned. First, we have to simplify Eq. (PT|); 
then apply similar arguments to which will 
allow to find 7. 

We are interested in the region of z = (zi + 
z 2 )/2 behind the barrier. We can make this re- 
gion rather small, compared to the region away 
from the barrier, by choosing potential U(z) to 
be infinite for z < zq. Here zq is an arbitrary, fi- 



rs 



nite point behind the barrier. This procedure 
does not restrict generality of this argument, 
since the actual region behind the Coulomb bar- 
rier of a fusing particle is much smaller than the 
region away from the barrier. Then, we can ap- 
proximate V(q±, Z1—Z4) by K(q_i_, Z—Z4), as well 
as V(q ± ,z 2 - z 3 ) by V(q±,z - z 3 ) in Eq. ([!!]). 
As we explained above, the main contribution 
to the integral over z 3 , z^ comes from the region 
of z > z* (the region away from the barrier). 



In this region, the functions g b (z 3 , k), gbiz^, k), 
entering G ( b ° )+ ~ , G ( b ° y+ in Eq. (0), can be well 
approximated by exp(ikz 3 ), expiikz^). After this 
approximation, the region of integration can be 
extended from —00 to 00, which will introduce 
only an error in the numerical coefficient on the 
order of unity. By using the explicit forms of 
G^~ + , G^ + ~ from Eq. (^), we can carry ex- 
plicitely integrations over frequencies. Then we 
find 



-4iri 



dp 1 <iq 

(2?r) 3 (2tt) 3 J (2tt) 



Vfflwp - e b pi + e^_ q - epx_ q x >fc ) 
In a similar manner, we find the answer for the 



(33) 



dpi <iq 

(2tt) 3 (2tt) 3 J (2tt) 



(34) 



The expression for 7 follows from Eqs. (|3^, |33], 
34|) . We see that its dependence on z is deter- 
mined by functions g(z,k), which are exponen- 
tially suppressed behind the barrier. This com- 
pletes the proof and allows us to neglect the "in- 
homogeneous" term in the Eq. (|26|). 

The solution of the "homogeneous" equation, 
G~f + , involves only one factor of W(k), since 
G 



7 

j + oc gg*, see Eq. (p6|) . Therefore, to obtain 
K we will find only the solution to the "homo- 
geneous" equation. We supplement this homo- 
geneous equation with the boundary condition, 
G h from Eq. ( p7|) imposed at z = z* away from 
the barrier. 

Since the largest contribution to K is made 
by Gj + one can see, that we can use Eq. (||) to 
find K if we substitute Gj + instead of G^~ + in 
(|). Therefore, we find Gj + with dependence on 
Zi, Z2 from the very begining by acting on Eq. 
with the operator 



g(z 1 ,k 1 )g*(z 2 ,k 2 



dk\ dk 2 
27T 2n 



(35) 



We use the boundary condition, Eq. (|27|) and 



obtain the following answer for Gj + : 

in F (uj) I g^z^b^uj - e^ k )gl(z 2 )dk (36) 



(37) 



where gu{z) = g(z, k) and we used for 8^ 



2ImG* , (k) 



(^-e q±fc ) 2 + (l^A k )f 



Comparing Eq. (|j) with (ffij), fl37|) we see, 
that the effect of collisions, which makes the 
largest contribution to K can be described as 



8{u 



<5 7 (u; 



e q ± fc) ; as can be expected 



on the intuitive grounds. Now we substitute 
Eqs. fl3"6D, ( |3"T| ) into (|5|) and integrate / du}/2n by 
using the result of H for large momenta, q 2> qr'- 



2^ 



n M (u)5^(u - e qxfc ) = n M (e^ k ) + 8n 7 (q±k), 



where <5n 7 (q_i_/c) is the "quantum tail". 

To describe the tunneling of particles mi and 
m 2 we perform the same steps as those leading 
to Eq. (]25|) . Note that this procedure implies 
quantization of the hypothetical particle with 
the mass equal to the reduced mass of fusing 



G 



particles. Strictly speaking, this approximation 
is correct in the limit when one particle is light 
and the other is heavy. Keeping this in mind, we 
arrive at the answer for K: 



show in Appendix A, the width, 7, can be ap- 



proximated as hu v 



Therefore, this means 



K — n 



duj p f°° 4irp 2 dp S(e t 



2tt 



(2*) 



p 

— W(p)n M (uj)S^(uj 



(3? 



We see that it involves integrals over both the 
energy and momentum. The width of the res- 
onance, 7(o> p , p), is, in general, a complicated 
function of energy and momentum. It obeys un- 
der certain assumptions the integral equation || 
which follows from the Dyson diagram expan- 
sion. More simplified calculations, based on the 
use of the zero-th order Green functions, lead to 
the analytical answer found in [Q (see Appendix 
A). We use this answer for 7 0|| and obtain 



K 



n 



00 47ip 2 dp S(e p ) 
(2vr) 3 e p 



[%(p) + Srijip)) = K M + K n 



(39) 



Finally, we evaluate S R through the collision fre- 
quency, v T , carry out the integral in (^) and 
form the ratio r 12 = K 1 /K M (see Appendices A 
and B): 



that in the limit of negligible width, the correc- 
tion to the rate becomes very small, and the re- 
action rate coincides with the usual Gamow rate. 
This is expected, since the smaller the width, the 
closer the particle to the mass shell. And the 
Gamow rate is derived for on-shell particles. On 
the other hand, the ratio becomes large, if hv p /T 
is not very small, yet still Tiu p /T < 1. It is this 
condition, which allows us to use the result of 
H for 7 in Eqs. (|59"D, ([II]). In the opposite case, 
hv p /T > 1, corresponding to more strongly cou- 
pled plasma, the theory presented here does not 
apply. 

Consider a DT plasma of density p = 10 g/cm 3 
and temperature T = 0.1 keV. The ratio 
TiVpjT = 0.0075 for these conditions. Then 
Eq. (|40|) will lead to the rate accelerated by 
five orders of magnitude as compared to the 
conventional answer. We find the ratio tbt = 
Kl T /K% T = 1.8 • 10 5 , with Kp T = 2.6 ■ 
10- 30 cm 3 /sec, K^ DT = 4.8 ■ 10- 25 cm 3 /sec. 

It is also interesting to consider how the reac- 
tion rates are changed under conditions relevant 
to astrophysics, e.g. in the Sun interior. Modi- 
fication of several reaction rates and their influ- 
ence on neutrino fluxes from the Sun was con- 
sidered in ||. Here we give examples of p + Be 7 
and p + p reactions, occuring at the center of 
the Sun, r = 0, where the temperature is T = 



ri 2 



3l9 ^ J2 huT ( mco11 ^ ( (40) 1.336 keV and the density p = 156 g/cm 



ac- 



'12 



2 72 j,„\ 1/3 



Here m co u = m r mj/{m r + m 

m 1 m 2 /(m 1 +m 2 ), r 12 = 3(vr/2) 2 / 3 ( 10 °^" 
Aw = AiA 2 /(Ai + A 2 ), Tfc is the temperature in 
keV; Ax, A 2 are the atomic numbers of tunnel- 
ing particles, rrij is the mass of the background 
particles colliding with tunneling particles. Sum- 
mation is over all species of colliding background 
particles, including tunneling species. This re- 
sult, (|39|) , (|40|) , supports the idea of averaging 



cording to the solar model of [21 
for 7l7 /T is %v T /T ~ 0.0056 



The estimate 
while the rates 



and the ratio, r 17 , are 



Kl 7 = 4.0 ■ 10 



-30 



9.6-10 



cm 3 /sec, ri7 = 4.17 



36 cm 3 /sec, 
10 5 . 



of av over momentum distribution postulated in 

Note an important feature of the result (flO]): 
this ratio goes to zero as hu p /T — > 0. As we 



In case of heavier elements and not very small 
hux/T, the ratio becomes much higher. On 
the other hand, the corresponding numbers for 
the p + p reaction under the same conditions 
are: j/T ~ twr/T ~ 0.00054, Kg = 1.4 ■ 
10~ 43 cm 3 /sec, r n = 3 ■ 10~ 3 . The reaction rate 
is essentially the same. Note that Tivt/T ~ 
0.00054 is rather small. This confirms the con- 
clusion of the discussion above, that small 7 lead 
to insignificant changes in the reaction rates. 
Note, however, that we need to have some gauge 



7 



to determine which 7 is small. Naturally, such 
a gauge should be a ratio of the height of the 
Coulomb barrier to the average particle energy. 
As we see from Eq. fl40D this ratio enters the an- 
swer for r 12 through r 12 . 

In conclusion, we derived the momentum dis- 
tribution in equilibrium through real-time tech- 
nique of Keldysh and Korenman which coin- 
cides with the distribution, obtained in imagi- 
nary time |4]]. This distribution takes into ac- 
count off-shell effects which lead to the quantum 
tail with power-like dependence on momentum 



in contrast to the distribution for on-shell par- 
ticles which is Fermi-Dirac or Maxwellian. We 
calculated the rates of fusion reactions for parti- 
cles in a weakly coupled plasma and showed that 
off-shell effects may significantly enhance the re- 
action rates. 
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APPENDIX A 

CALCULATION OF THE MASS OPERATOR, AND THE DISTRIBUTION 

FUNCTION OVER MOMENTUM 



Galitski and Yakimets found how off-shell pro- 
cesses change the distribution function over mo- 
mentum in equilibrium, using imaginary time 
technique |4[]. In this Appendix we show how 
to find the same result through real time dia- 
gram method. We also generalize it slightly by 
allowing for collisions between different types of 
particles. 

Dyson equations for the function G h can be 
written as 

R,t;w,p) = 

(41) 

To simplify notation we do not write explicitely 
the variables the functions £ and G depend 
upon. We set to zero the LHS of this equation 
in equilibrium and obtain 

G-+ = ^-(G R -G A + G- + ) , (42) 

which follows from the relation between Green 
functions [13|]. The result for G~ + follows by 
substituting Eq. (g|) into fl2|): 

ImL R 




G' 



—2mnp(uj)- 



[oj-e p - ReZ R ) 2 + (ImZ R ) 2 

(43) 

We use this equation to find the distribution 
function over momentum which is 



/(P) 



~27 



G~ 



(44) 



First we calculate the mass-operator H R . It can 
be expressed using the relation 



2 V 



2 V 



(45) 



The first term in this equation corresponds to 
the real part ReH R , while the second term is 
imaginary part ImH R . Since it is the imaginary 
part of Tj R which is responsible for the quantum 
tail 0] we calculate only Im'E R and ignore ReL R 
altogether. 

The first two terms of the series for the oper- 
ators S a/3 are shown in the Fig. [I]. One has to 
assign all possible combination of signs — , + to 
the vertices, keeping in mind that the sign does 
not change along the dashed line M. 



p.- Q 



P- Q 



P- Q 



P- Q 




P+Q 



a) b) c) 

FIG. 1. First and second-order terms in the diagram ex- 
pansion of T, R . 



8 



In a more rigorous approach, based on the 
resummation of the series of graphs, one can 
show flTDU9U?H that, solid lines in the diagrams 
of Fig. [i] correspond to exact Green-functions, 
while the dashed lines represent scattering am- 
plitudes of the particles in the media. However, 
such resummations lead to intractable non-linear 
integro-differential equations. 

As a first approximation we assume that, the 
solid lines correspond to zero-order Green func- 
tions, while the dashed lines correspond to the 
interaction potential between the particles. 

These diagrams represent quite different phys- 
ical phenomena. Diagram "a)" describes the 
screening effects and formation of the quasipar- 
ticles in the plasma. Therefore, an imaginary 
part of this graph represents decay of the quasi- 
particles. This term is conventionally put into 
the LHS of the kinetic equation, signifying that it 
has nothing to do with the collisional integral J7| . 
Diagram "b)", on the other hand, if evaluated 
with zero-order Green-functions, gives the usual 
Boltzman collision terms fl3lj7l1. Therefore, the 
imaginary parts of the diagrams "b)" and "c)" 
describe damping of the oscillators used for the 
second quantization of the interacting particles. 

We evaluate the "b)" and "c)"-type of dia- 
grams for assuming that collisions occur 
between particles of type a and b. Therefore, we 
assign the index a to the Green function, rep- 
resented by the horizontal line in the diagram 
"b)", and the index b to the Green functions 
of the loop. Note that since the diagram "c)" 
corresponds to the exchange processes only sin- 
gle index can be assigned to its Green functions. 
This obviously reflects the fact that only iden- 
tical particles are subject to exchange. Taking 
into account these remarks, we write the analyt- 
ical expression for the diagramm b) in the form 

- iE-+ = J G-+(P - Q)G+-(P 1 )G; + (P 1 - Q) 

I^WW' (46) 



where capital letters denote four-momentum, 
e.g. P = (lu p , p) and V(q) is the interaction po- 
tential between the particles. We substitute here 
equilibrium zero-th order Green functions [[HJ 

G~ + (lu p p) = 2ni n p 5(tu p - e p + /i), (47) 

G + ~(u p p) = — 2-7TZ (1 - n p )5(uj - e p + //). (48) 

Analogously, we obtain an expression for the di- 
agram c). 

For a moment, we omit the indices a and b to 
simplify the formulas, but write them explicitely 
in the final formulas. Then the result for E h is 



E-+K,p) = (-2™)/^L/ 
(2V q 2 - V q K 



dpi f <iq 

(27p J (2^)3 



q-p-pi I ^p-q^pi-q 



<K^p - e P -q + e Pl -e P i-c 



(49) 



Proceeding in a similar fashion with E + and 
using Eq. (f45|) we obtain the final result for 

f dq dpi 

1 (uj p ,p)=mJ _ _F q (2y q - F q _ p _ pi ) 

[n p _ q (n Pl _ q - n Pl )+ n Pl (l - n Pl _ q )] 
5(^ p -e p _ q + e pi -e Pl _ q ) (50) 

This answer coincides with the result of p2|fl 
obtained in the imaginary-time technique. 

Note that we could improve this relation by 
substituting <5 7 of Eq. (|37| ) instead of the delta- 
function. This, however, would lead to an inte- 
gral equation for 7, and its first iteration will be 
just Eq. (0). 

Now we can find the distribution function over 
momentum from Eqs. (0), (0), P3| ) 



/(P) 



d(jJr 



n F (u)~f(u p ,p) 



(27r) 2 K-e p ) 2 + 7 2 K,p) 



(51) 



by substituting Eq. (0) into Eq. © *. The 
integral / dou can be rewritten as a sum of two 
terms 



*From now on we follow closely the work of | 
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£ ° d 



-oo 
oo 



2n 

d(J v 
~2t7 



(52) 



CO 



where eo is the characteristic energy. In the de- 
generate case, T <C we have eo = /i, in the 
nondegenerate case, T ^> fi and e = T. We 
define p as e = pi /2m. 
Consider the case of large momenta, p ^> p . 



The function <5 7 (u; p 



e p ) has a maximum at 



lu* = e p ^> e , so that co>* G [e , oo], and the con- 
tribution of the second integral can be approxi- 
mated as rip. In the region of integration of the 
first integral in Eq. (|52|) we can set np(uj) ~ 1. 
Then Eq. ([52] ) becomes for p 3> po 



/(P) 



cicup 7(^ P ,P) 
2vr (u p - e p ) 2 



n p + 5n 7 (p), 
(53) 



We carry out the integral in the second term of 
(|53|), use Eq. (|50"D and obtain 



5n 7 (p) 



(iq dpi 



(2tt) 3 (2tt) 3 



Vq(2Vq - K 



q-p-pi, 



[np_ q (n pi _ q -n pi ) +n Pl (l 



n 



pi-q 



)] 



e p-q 



pi 



e pi-qJ 



with the condition on momenta 



oo < e p q 



-Pi 



31 -q < e 



(54) 



(55) 



due to the integration of the delta-function from 
7(w P ,p). 

The integral in Eq. (|54]) should be carried out 
numerically. We can, however, obtain a reason- 
able estimate in the case of large momentum, 
p ^> po- Observe that the first two terms in the 
brackets [..] should have the same dependence on 
p. To obtain an estimate consider the product of 
the first term in (..), V 2 , and the first term in [..] 
which involves the product n p _ q n pi _ q . We de- 
note it ii(p). Note that there is a region of non- 
exponential contribution to the integral, which 
is limited to the region of momenta 

|p - q| ~ |pi - q| ~ p ,p ~ p 1 ~ q > p (56) 



The inequality (|5^) is satisfied which, in turn, 
means that the dominant contribution to this 
product can be approximated as 



(a) (6) 
2 n p -q n pi-q 

el 



, (57) 



where e p = p 2 /2m a h and we reinstalled indices a 
and b. Consider non-degenerate case. 

We change variables, p — q = q', pi — q = p[ 
and obtain an estimate 



2 2 
e a e b~ 



(2 m ab f 

p8 



\ 3/2 / x 3/2 



2(3tx 



2(3n 



We substitute the relation 




2/37T, 



3/2 



(58) 



(59) 



(60) 



^(eAl^n^ (61) 



Since the physical reason for non-zero j(uj p , p) is 
collisions, we rewrite (|H|) in terms of the colli- 
sion frequency 



as 



7TX3/2 

2) 



n^(e a e b ) 2 m:^'T 



1/2^-3/2 (2m ab T) 3 / 2 



p° 



2 hi/pT pt n 



(b) 



71 



3/2 72 



73> ' 



p p- T 



(62) 



(63) 



where px = {2m ab T) 1 ^ 2 . This is essentially the 
result of [|] for the degenerate case, sligtly gen- 
eralized to the case of colliding particles with 
different masses. 

In the non-degenerate case we can keep only 
the term linear in n p in Eq. ([54]) , because the 
quadratic terms are smaller. We denote this 
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term by I3. Since there is no exponential sup- 
pression in the g-integral for this term, we have 
to evaluate the limits of the g-integration care- 
fully. We can obtain these limits from Eq. (j55|). 
After some algebra it leads to 



2m ab 



— qx 



> e n — T, 



(64) 



_E_ + Pi 

m a m b 

where x is the cosine of the angle between q and 
— + — . We resolve this inequality in the limit 
of large p 3> p?, neglecting all terms on the or- 
der of p\ in comparison to p 2 . Note that since 
the term in quesiton, J 3 , is proportional to n pi , 
the contribution from the region pi > p? is ex- 
ponentially suppressed. Since we are interested 
in the non-exponential contribution only, we can 
assume that p\ ~ pt in Eq. (|64]) . Then we finally 
arrive at 



m ab 

q> q* = p{l 



'1 + 



m ab 



(65) 



To obtain an estimate for J 3 , we carry out the 
integral of the linear term in n Pl from Eq. ([54]), 
multiplied by the first term in (..), V3: 



1 



*s(p) 



(e a e b 



2tt 



dpi J dx 

-1 q* 



q 2 dq 



qx 



p 1 pi 



(66) 



After performing the integral over x, we are left 
with the integral over y, which is the cosine of 
the angle between p and pi. The y- dependent 
term becomes 



/ 



dy 



1 



V 2m at 



q 



^ + 

mi 



+ 



m a m b 



(67) 



We carry out this integral, neglect all the terms 
of the order of p\ in comparison to p 2 and obtain 



2m 2 /(q p ). The remaining integrals over / dpi 



00 

and / dq are trivial. Using Eq. (|65| ) we get the 

q* 

result for J 3 



m lf m a\ 3 (e a e fe ) 2 n (b) e /3/ia 
3(2tt) 3 ^Uw 



1 + J1 + ^ 

m ab 



Now we can see that different terms in Eq. ( p4|) 
have a different dependence on p. Note however 
that the answers for these terms are really es- 
timates. To find the correct answer one should 
solve non-linear equation for 7 []5] , which is pos- 
sible only through numerical methods. 

It is also possible to obtain an estimate 
through dimensional analysis which leads to p[ 



7 



2ImY, R ~ 2hi 



(6* 



This expression has to be substituted into 
Eq. (|53|). Note that we have to carry out the 
integral over energy from — 00 to T. Since we do 
not know how 7 behaves at negative energies, we 
can only make an assumption that, this region 
does not contribute significantly to the integral, 
and the main contribution comes from u p ~ T. 
Then the integral is equal roughly to 



I{p) 



hv p T n (a) 



We can simplify it to 
J(p) 



hu p T r ^ a 



(69) 



(70) 



which is the result for the non-degenerate case 
found in 0. We will use this estimate in our 
calculations of the reaction rate because of its 
simplicity. 



APPENDIX B 

REACTION RATE: CALCULATION OF THE INTEGRAL OVER 5n 7 (p) 
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In this Appendix we find the rate of 
Eq. (p9|). Observe, that in the main text we 
used approximation of a hypothetical particle 
with the mass equal to the reduced mass of fus- 
ing particles, m r . It is this particle which under- 
goes collisions with the background, leading to 
its power-like distribution over momentum. We 
will need the result fl69|) for these particles. We 
write as 



Hp) 



2-Kil \m coU ) 



taking into account that 



Pt 



\m co uJ 



(71) 



(72) 



Therefore, we substitute Eq. (|TT|) into Eq. ([39] 



and obtain K~» 



( ^-Y /2 (2m r y hMmc ° ll)T 
\m co uJ 



2tt 



S{t*){2m coU T) 



3/2 



-wp G /p 



—dp, 



(73) 



where S(e*) is the astrophysical factor at some 
fixed value of energy. The integral is trivial and 
is equal to 



-Ttpa/p 



1 



■ dp 



24 

7T 5 P 5 G 



(74) 



For the sake of completeness, we also find here 
the usual Gamow rate: 



K 



M 



4ttS{6*) 



-TTp G /p-p 2 /2m r T 



pdp. (75) 



The result, obtained by using the steepest de- 
scend method, is 



K 



M 



S(e* 



m r T)r 



l/2 p -r 



(76) 



Strictly speaking, the astrophysical factors in 
Eqs. ([T3|) and ( |76| ) need not be the same. In fact, 
one should consider a generalized astrophysical 
factor, S(ujp, e p ), depending both on energy and 



momenta. Such factors would take into account 
the influence of the off-shell effects on nuclear 
transformations. However, as far as we know, 
they have not been studied. Therefore, we as- 
sume that the factors in Eqs. fl73|), flT6| ) are the 
same. 

With this in mind, we can form the ratio 



[2. 
[3 
[4 

[5 

[6 

\-\ 
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[11 

[12 
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[is; 
[17; 
[is; 
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' 2 2 tt 3 /2 t \m r J 
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